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ABSTRACT 

We present evidence for Quasi-Periodic Oscillations (QPOs) in the recurrent outburst emission 
from the soft gamma repeater SGR 1806-20 using NASA's Rossi X-ray Timing Explorer (RXTE) 
observations. By searching a sample of 30 bursts for timing signals at the frequencies of the QPOs 
discovered in the 2004 December 27 giant flare from the source, we find three QPOs at 84, 103, and 
648 Hz in three different bursts. The first two QPOs lie within ~ 1 a from the 92 Hz QPO detected 
in the giant flare. The third QPO lie within ~ 9 a from the 625 Hz QPO also detected in the same 
flare. The detected QPOs are found in bursts with different durations, morphologies, and brightness, 
and are vindicated by Monte Carlo simulations, which set a lower limit confidence interval > 4.3cr. 
We also find evidence for candidate QPOs at higher frequencies in other bursts with lower statistical 
significance. The fact that we can find evidence for QPOs in the recurrent bursts at frequencies 
relatively close to those found in the giant flare is intriguing and can offer insight about the origin 
of the oscillations. We confront our finding against the available theoretical models and discuss the 
connection between the QPOs we report and those detected in the giant flares. The implications to 
the neutron star properties are also discussed. 

Subject headings: pulsars: individual (SGR 1806-20) — stars: flare — stars: neutron — stars: 
oscillations — stars: magnetic field 



1. INTRODUCTION 

Soft Gamma Repeaters (SGRs) and Anomalous X-ray 
Pulsars (AXPs) are widely believed to be magnetars, 
i.e. neutron stars with inte nse magnetic field of the or- 
der of 10 14 - 10 15 Gauss ([Duncan fe Thompson! [1991 : 
IWoods fe Thompsonll2006| ) . Magnetars are characterized 
by the repetitive emission of bright bursts of low-energy 
gamma-rays with short durations (~ 0.1 s) and peak lu- 
minosities reaching up to 10 41 erg s~ x , well above the 
classical Eddington luminosity L e dd ~ 2xl0 38 erg s" 1 
for a 1.4 Mq neutron star. Magnetars are also persis- 
tent X-ray pulsars with slow spin periods (P ~ 2 — 12 
s) and rapid spin-down rate (P ~ 1 ms/yr). Less fre- 
quently, magnetars produce enormously energetic giant 
flares with peak luminosities of the order of 10 44 - 10 46 
erg s _1 , followed by pulsating tail lasting for hundreds 
of seconds. At pres ent, there a r e 18 m agnetars: 8 SGRs 
and 10 AXPs0 (see iMereghettil (|2008f ) for a review). 

The magnctar model postulates that the dominant 
source of energy of magnetars is their intense magnetic 
field. The persistent X-ray emission is attributed to mag- 
netospheric currents driven by twists in the magnetic 
field. The short bursts are produced as a result of out- 
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ward propagation of Alfven waves through the magneto- 
sphere. The giant flares are triggered by sudden catas- 
trophic reconfiguratio n of the magnetic field that causes a 
large scale crustquake dThompson fc Duncanlll995l 119961 . 
[20011 ISchwartz et al.ll2005H . 

On 2004 December 27, SGR 1806-20 emitted its first 
recorded giant flare that was observed by a number of 
satellites as the most powerful giant flare from a mag- 
netar. The flare reached a peak lu minosity of the or- 
der of 10 46 erg s~i in the first 0.2 s (| Hurley et al.l[2005t 
iPalmer et a.l.ll2005h . causing a strong perturbation in the 
Earth's ionosphere and satur ation of most satellite de- 
tectors (jTerasawa et al.|[2005f ). The initial spike was fol- 
lowed by ~ 400 s tail modulated at the 7.56 s spin pe- 
riod of the SGR. Using R XTE Propor t ional Counter Ar- 
ray (PCA) observations, llsrael et al.l (|2005D discovered 
three QPOs at 18, 30, and 92 Hz in the decaying tail. 
With Ram aty High Energy Solar Spe ctroscopic Imager 
fRHESSn. lStrohmaver fc Watts! ([2005h confirmed the 18 
and 92 Hz QPOs and report ed two additional QPOs at 26 
and 626 Hz. Subsequently, I Watts fc Strohmaverl (|2006Q 
re-analyzed the RXTE observations and confirmed the 26 
and 626 Hz QPOs detected with RHESSI and reported 
additional QPOs at 150 and 1840 Hz. Two more QPOs 
a t 720 and 2384 Hz wer e repor ted at a lower significance. 

iStrohmaver fc Watts! p005T ) also studied the RXTE 
observations of the 1998 August 27 giant flare from 
SGR 1900+14 and discovered four QPOs at 28, 53, 84, 
and 155 Hz in the flare's pulsating tail. In both sources 
the QPOs were detected in part of the tail and were 
only seen at a specific rotational phase. They were 
identified as a sequence of toroidal modes associated 
with seismic vibrations in the neutro n star crust. Other 
model s were also s u ggested, e.g. iGlampedakis et al.l 
(|200l) : ILevinl (|200l : iBeloborodov fc Thompson! (|2007fl ; 
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Samuelsson fc Anderssonl (120071) ; iMa et all (|2008t) ; 
Shaisultanov fc Eichlerl (|2009l) . 



Motivated by these findings, we undertook the task of 
searching for similar signatures in the typical and more 
common emission of magnctars: the recurrent bursts. 
Here we present the results from SGR 1806-20. Our re- 
sults from SGR 1900+14 are forthcoming. 

2. DATA ANALYSIS & RESULTS 

RXTE first observed SGR 1806-20 during 1996 Novem- 
ber 5-18. We obtained the PCA observations in 2-60 
keV and reduced the data with HEAsoft 6.6 to obtain 
the burst light curves. The observations with the most 
bursts were taken in the event-mode configuration with 
time resolution At ~ 125/is and 256 energy channels. 
All 5 proportional counter units (PCUs) of the PCA 
were fully functional. Photons from all PCUs and all 
detector anodes were combined to maximize the signal- 
to-noise ratio and the artificial event-mode time mark- 
ers were removed. We screened out very faint bursts 
that do not have sufficient counts and very bright bursts 
that would cause significant pile-up and dead-time ef- 
fects and selected 30 bursts for our timing analysis. We 
binned each burst at the smallest time resolution al- 
lowed by the data configuration (125 /is) and applied 
Fast Fourier Transformation (FFT) to the binned pho- 
ton arrival times to obtain a power spectrum for each 
burst. The frequency resolution A.P = ^txN b - — ranged 
from 1.72 to 14.63 Hz, depending on the burst duration, 
and the Nyquist frequency was v n yquist = jki ~ 4096 Hz. 
The power spect ra were normalized using Leahy normal - 
ization method (iVan der Klislll989l: Hsmel fc Stellalll996D 
and we searched each burst power spectrum for tim- 
ing signals in the frequency range of the QPOs found 
in the 2004 giant flare. Candidate QPOs were fitted 
to a Lorcntzian function to calculate their properties, 
and their significance was investigated using Monte Carlo 
simulations. 

In a short (~ 0.09 s), moderately bright burst (total 
counts ~ 1.13 x 10 3 ), we find a QPO candidate around 
84 Hz with single trial null probability of 3.4xl0 -5 (see 
Fig. 1 top panel). The QPO has a centroid frequency 
v Q = 84.31 ± 0.66 Hz and a half width at half maximum 
a = 8.69 ± 0.78 Hz, corresponding to a coherence value 
Q = v /2a of 4.84. 

In another short, bright burst with a gradual rise (du- 
ration ~ 0.09 s and total counts ~ 3.07 xlO 3 , see Fig. 1 
middle panel), we detect a candidate QPO at 103 Hz at a 
null probability of 1.43xl0~ 6 . This QPO has a centroid 
frequency v a = 103.44 ±0.42 Hz, width a = 10.65 ±0.42 
Hz, and coherence value Q ~ 5. 

In the bottom panel of Fig. 1, we show a long (> 0.1 
s), bright burst (total counts ~ 1.37 xlO 4 ) with a pair 
of QPO candidates at around 648 and 1096 Hz with a 
single trial probability of 2.33xl0~ 7 and 1.37xl0~ 5 , re- 
spectively. Fitting each to a Lorentzian function yields 
v Q = 648.5 ±0.15 Hz, a = 2.5 ±0.09 Hz, and Q = 130 for 
the first QPO and v Q = 1095.88±0.21 Hz, a = 2.88±0.17 
Hz, and Q = 190 for the second. 

In the other bursts we did not find significant QPOs 
at the previously detected frequencies. Extending our 
search to other frequencies, we detect two high frequency 
QPO candidates in the short, relatively dim burst (du- 
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Fig. 1. — Light curves, power spectra, and simulated power spec- 
tra of the three bursts that showed evidence for QPOs. Under each 
burst is its power spectrum (left) and the accumulation of 50,000 
power spectra (right) obtained with Monte Carlo simulation. The 
top burst shows a QPO at 84 Hz (Ai> = 11.38 Hz), the middle 
burst shows a QPO at 103 Hz (Au = 11.38 Hz), and the bottom 
burst shows a pair of QPOs at 648 and 1096 Hz (Au = 4.27 Hz). 
The dashed horizontal line represents the baseline of the single 
trial 95% confidence level. The simulated power spectra reveal no 
significant structure at or in the vicinity of the detected QPOs. 
The plotted light curves are binned at 2 ms. The bursts occurred 
on 1996 Nov. 18 at 06:18:44.446 UTC (top), 08:57:02.083 UTC 
(middle), and 04:41:53.480 UTC (bottom). 



QPOs in the Reccurent Emission from SGR 1806-20 



3 



ration ~ 0.07 s and total counts ~ 700) shown in Fig. 2 
(upper burst) at 1229.8 ± 0.80 and 3690.74 ± 0.67 Hz. 
Forming the null hypothesis results in a single trial prob- 
ability of 3.4xl0 -5 for each feature. We estimate that 
the 1230 and 3690 Hz features have a = 8.00 ± 0.48 and 
8.37 ± 0.79 Hz respectively and corresponding coherence 
values of 76 and 220. Another high frequency QPO can- 
didate around 2785 Hz is found in the long, intense burst 
(duration ~ 0.63 s and total counts ~ 2.00 x 10 4 ) shown 
in Fig. 2 (bottom burst), at a single trial probability of 
7.96xl0" 6 . The QPO is centered at v Q = 2785.36 ±0.20 
Hz and has a width a = 0.75 ±0.19 Hz and coherence 
Q = 1860. 

Since our strategy was to search for QPOs at the pre- 
viously reported frequencies, the number of trials for the 
84, 103, and 648 Hz QPOs takes into account the num- 
ber of bursts (30 bursts) and the number of frequency 
bins between the QPOs we detect and those discovered 
in the giant flare (0.72 bins for the 84 Hz QPO, 0.96 
bins for the 103 Hz QPO, and 5.51 bins for the 648 Hz 
QPO). This brings their null probabilities to 7.34xl0~ 4 , 
4.14xl0~ 5 , and 3.85xl0~ 5 , respectively. For the high 
frequency QPOs at 1096, 1230, 2785 and 3690 Hz, we 
had no prior knowledge of these frequencies and we must 
account for the total number of frequency bins in the 
power spectra as well as the number of bursts. The num- 
ber of frequency bins accounted for in the case of the 
1096, 1230, 2785 and 3690 Hz QPOs is 960, 280, 2376 
and 280 bins respectively. This brings the chance prob- 
ability to higher than 10 -3 , thus ruling out these QPOs 
on the basis of the null hypothesis alone. 

A more cogent approach to assessing the statistical 
significance is through utilizing Monte Carlo simulation. 
For each of the above bursts, we generated a simulated 
light curve at the same time resolution (125 /is) by fit- 
ting the burst light curve to a high-order polynomial 
function. We then determined the number of photons 
in each time bin by seeding a Poisson distributed ran- 
dom number generator whose mean value is the number 
of counts we obtain from the fitting function. We used 
this approach to obtain 50,000 simulated light curves for 
each real burst and we generated a power spectrum for 
each simulated burst using the same technique we used 
to find the QPOs in the bursts above. In Fig. 1 and 2, 
the right panel under each burst shows the accumulated 
power spectra of the 50,000 simulated burstfl We find 
none to have shown a structure at or near the frequen- 
cies of the candidate QPOs with similar or higher power, 
yielding a null probability lower limit of < 2 x 10~ 5 (or a 
confidence level interval > 4.3 a) for the QPOs reported 
above. 

The statistical significance arrived at from the two in- 
dependent approaches above gives a secure detection for 
the 84, 103, and 648 Hz QPOs and a marginal signifi- 
cance for the kHz QPOs. 

3. DISCUSSION 

Our search for QPOs in the recurrent bursts from 
SGR 1806-20 revealed QPOs of interesting similari- 
ties to those discovered in the gia nt flares from the 
source and from SGR 1900±14 (jlsrael et al.l [20051: 

2 The number of simulated bursts was set by the available CPU 
and memory resources to carry out the Monte-Carlo simulation. 
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Fig. 2. — Light curves, power spectra, and simulated power spec- 
tra of the two bursts that showed evidence for QPOs at high fre- 
quencies. Under each burst is its power spectrum (left) and the 
accumulation of 50,000 power spectra (right) obtained with Monte 
Carlo simulation. The top burst shows two QPO at ~ 1230 and 
~ 3690 Hz (Au = 14.63 Hz). The bottom burst shows a QPO at 
~ 2785 Hz (Ai> = 1.72 Hz). The dashed horizontal line represents 
the baseline of the single trial 95% confidence level. The simulated 
power spectra reveal no significant structure at or in the vicinity 
of the detected QPOs. The plotted light curves are binned at 1 ms 
(top) and 2 ms (bottom). The bursts occurred on 1996 Nov. 18 at 
06:18:44.134 UTC (top) and 12:18:01.487 UTC (bottom). 



iStrohmaver fc Watts! [20051: iWatts fc Strohmaverl [2001 . 
We found three intriguing QPOs at 84.31 ±0.66, 103.44± 
0.42, and 648.5 ± 0.15 Hz. The 84 Hz QPO lies within 
~ 1 a from the 92.5 ± 0.2 Hz QPO from the same source 
and coincides with the 84 Hz QPO discovered in the SGR 
1900±14 giant flare. The 103 Hz feature lies within ~ 
1 a from the 92.5 ± 0.2 Hz feature from the source. The 
high frequency QPO at 648 ± 0.15 Hz is within ~ 9 a 
from the 625.5 ±0.15 and 626.5 ±0.02 Hz QPOs detected 
from the source with RXTE and RHESSI, respectively. 
Those QPOs are detected in bursts with different dura- 
tion, luminosity, rise time, and light curve morphology, 
and thus cannot be due to artifacts in the light curves. 
Their statistical significance is further confirmed by the 
Monte-Carlo simulations. 
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Torsional shear mode mechanism on the neutron star 
crust is a promising framework for interpreting magnetar 
QPOs since the observed low frequencies fit comfortably 
within theoretical models predictions and are associated 
with a sequence of I modes, while the high frequency 
QPOs are interpreted as modes with higher order ra- 
dial displacement eigenfunction (n > 1). As a conse- 
quence of the twisting of the magnetic held lines and the 
outward propagation of the Alfven waves, stresses build 
up in the neutron star crust, which is very likely to ex- 
cite toroidal vibrational modes in the neutron star crust 
(lRudermanlll991l; IRuderman et allll998t ISchwartz et al.l 
I2005Q . The excited frequency modes depend on the 
neutron star mass and radius, crustal composition and 
the magnetic field (lDuncanlll998l; lHansen fe Cioljl | 198tt 
McDermott et al.l 119881 iMessios et alJl200lH Pirofl200a 
Schumaker fe Thorndll983HStrohmaverlll991h . 



In t he context of global seismic vibrations, iDuncanl 
(|1998|) estimated the eigenfrequency of the fundamental 
toroidal mode 1=2, n=0 denoted 2^0 to be 



v{ 2 t ) = 29.8 



1.71 - 0.7M 1A R W 2 



0.87-Rio + 0.13Mi.4-R 



Hz. 



(1) 



10 



Where R w = R/10 km and M 1A = M/1.4M©. If 
a twisted magnetic field is embedded throughout the 
neutron star crust, one needs to account for the mag- 
netic field correction and the eigenfrequencies are modi- 



fied to v = 



1+ #- 



nl/2 



Hz, where v Q is the non- 



magnetic eigenfrequency, = \JAm\x ~ 4xl0 15 p\f 
Gauss and pu ~1 is the density in the deep crust in 
units of 10 g cm -3 . The higher order mode eigen- 
frequencies (n=0 & I >2) were shown to be scaled as 
v{it ) oc Ujl + l)] 1 / 2 (lHansen & Cioffi 1980) and were de- 
rived bv lStrohmaver fe Wattsl (|2005t) under the assump- 
tion that the magnetic tension boosts the field isotropi- 
cally. One must note that the degree to which the mag- 
netic field modifies the eigenfrequency is highly depen- 
dent on the magnetic field configuration and other non- 
isotrop ic effects which coul d alter this correction signifi- 
cantly (|Messios et al.|[200ll ). The eigenfrequencies of the 
higher order toroidal modes arc given by 



v{it a ) = v{gt ) 



1(1 + 1) 
6 



1/2 



1 1/2 



Hz (2) 



The QPOs at 92.5 and 30 Hz discovered in SGR 1806- 
20 giant flare are attributed to the v(^/t ) and 1^(2 i Q ) 
modes, respectively while the 626.5 Hz QPO was mod- 
eled as n=l toroidal mode where the choice of I is ir- 
relevant. In this context, the 84 and 103 Hz QPOs that 
we report correspond to v{et ) and v{st ) respectively. 
It was recently suggested that a number of peaks can 
be seen in the in the 78105 Hz frequency range of the 
power spectrum of SGR 1806-20 giant flare that could 
be attributed to magnetic splitting under the perturba- 

tive condition < 1 (|Shaisultanov fe Eichlerl I2009D . 

In that view, the two QPOs at 84 and 103 Hz may be 
the result of magnetic splitting of a degenerate mode. 
The 648 Hz QPO we report may be interpreted in 



terms of n=\ toroidal mode irrelevant of the choice of I. 
For such high frequency QPOs, the n >1 torsional modes 
in the neutron star crust is a more plausible interpreta- 
tion and they fall within the frequency ran ge suggeste d 
for the search for higher frequency QPOs bv lPirol (|2005| ). 

The detection of QPOs in magnetars has triggered a 
boost in the theoretical modeling that included attempts 
to construct more refined models that consider crust/core 
coupling, magnetic field geo metry, elastic prope r ties o f 
the crust and other effects. IGlampc dakis e t al.l (|2006D 
used a simple global MHD toy model to calculate the 
global magneto-elastic modes taking into consideration 
the magnetic coupling between the elastic crust and the 
fluid core under the assumption that the magnetic field 
is uniform. The model provides intriguing results and 
indicates that the system has a rich spectrum of global 
modes and it will naturally excite modes that minimize 
the energy transfer between the crust and the core. The 
excited modes frequencies are similar to those of the ob- 
served QPOs, where the modes that survive are those 
for which the coupling is minimum and have very weak 
amplitudes in the core. This model further supports the 
global seismic vibrations interpretation of the observed 
QPOs and it also accommodates for the 18 and 26 Hz 
QPOs seen in the SGR 1806-20 giant flare. A gen- 
eral relativistic formulation using the Cowling approx- 
imation, which does not include any correction for the 
magnetic field yielded a set of analytic estimates of the 
eigenfrequency modes based on numerical calculations 
(jSamuelsson fc Anderssonl 120071) . This model yields a 
different scaling with I than the previous [Z(Z-l-l)] 1 / 2 scal- 
ing, which leads to different mode assignments. Futhcr- 
more, the assumption of an isotropic magnetic field is 
inadequate as recent MHD simulations of the evolution 
of the magnetic field and its geometry in magnetars 
have indicated that magnetars have a strong dipole field 
component and a possible s trong internal poloidal field 
(jBraithwaite fc Spruitll2006f ). Since the QPOs we report 
and those observed in the giant flare are about 8 years 
apart, the observed frequencies can offer insight on the 
magnetic field decay of the neutron star. If we assume 
that the 103 and 92 Hz QPOs are attributed to the 1=7 
mode, we infer that the magnetic field strength of SGR 
1806-20 has decayed from 2xl0 15 in 1996 to 3.5 xlO 14 G 
in 2004. 

Magnetospheric oscillation models were suggested as a 
possible i n terpre tation for the origin of magnetar QPOs. 
iMa et all (120081 ) argued that the low frequency range 
QPOs could be modeled using the standing sausage 
mode oscillations of flux tubes in the magnetar corona. 
It is assumed that part of the plasma ejected during 
the giant flare is trapped by the magnetic field and 
then forms magnetic flux tube structures similar to 
what is seen in the solar corona. The tube oscillation 
model cannot accommodate the high frequency QPOs 
in SGR 1806-20. For the very high frequency QPOs, 
iBeloborodov fc Thompso'nl (|2007| ) suggested that the in- 
ner electric current in magnetar corona will induce very 
high frequency QPOs ~ 10 kHz where the corona self 
organizes itself into quasi steady state. Such high fre- 
quency QPOs remain to be detected while the lower fre- 
quency oscillations would involve the outer corona where 
the emission in this region is not intense enough to excite 
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the observed QPOs. The QPOs we report are detected in 
short (< 0.1 s), relatively dim bursts and in long (> 0.1 
s), bright bursts that would produce different signatures 
in the magnetosphcre. Since the short, weak bursts may 
not produce a fireball at all or be of magnetospheric ori- 
gin, this casts doubts on magnetospheric models as an 
origin for the QPOs in magnetars. 

A number of theoretical models postulated that AXPs 
and magnetars are neutron stars surrounded by fossil 
disks that were acquired during super nova collapse or 
during a common-envelope interacti o n (IChatteriee et al 
200tt IChatteriee fc Hernquistl I2000t IPerna et all 1200 " 



Wang et all (|2006ft reported evidence for a debris disk 



around the AXP 4U 0142+61, which led to the sugges- 
tion that the magnetar QPOs could be due to a possible 
disk as in the case of accreting neutron stars, which ex- 
hibit kHz QP Os that are believe d to originate in the ac- 
cretion disk (jVan der Klisl 120061 ) . The reported QPOs 
are detected during the bright X-ray emission of the 
bursts and are not seen in the decaying phase, thus ar- 
guing against arising from a possible passive debris disk 
around the magnetar. Furthermore, this interpretation 
is mainly ruled out by the fact that we have no evidence 
for the existence of a fallba ck debris disk around SGRs. 
iWatts fc Strohmaverl (|2007l) pointed out that for neutron 
stars that exhibit extremely high frequency QPOs, the 
inner disk radius of the debris disk is comparable to that 
of the neutron star. In the case of magnetars, the inner 
disk radius would have to be several solar radii to induce 
the observed magnetar QPOs, which rules out the modes 



of a debris disk mechanism as a viable interpretation for 
the magnetar QPOs. 

Levin (2006) argued against the global seismic vibra- 
tions as a viable interpretation for the QPOs in magne- 
tars and suggested that the QPOs could be associated 
with MHD elastic modes of the neutron star. However, 
the toy model employed by Levin assumed a rather sim- 
plified model for the geo metry of the neutron s t ar and the 
boundary conditions. IWatts &: Strohmaverl (|2007l ) of- 
fered a detailed critique and casted doubts about Levin's 
conclusions except for the crust /core coupling effect 
which needs to be well accounted for in future models. 

Observing QPOs in both the giant flares and the 
recurrent burst emission from magnetars is an impor- 
tant development in our understanding of neutron star 
physics. Further observations, analysis of archival data, 
and model developments are key to shedding further in- 
sight on the underlying phenomena. 
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